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Several mathematical views of phase-locking are developed. The classical Huyghens ap- 
proach is generalized to include all harmonic and subharmonic resonances and is found to 
be connected to 1 // noise and prime number theory. Two types of quantum phase-locking 
operators are defined, one acting on the rational numbers, the other on the elements of 
a Galois field. In both cases we analyse in detail the phase properties and find them 
related respectively to the Riemann zeta function and to incomplete Gauss sums. 
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1. Introduction 

In the crude sense phase-locking occurs whenever the erratic behavior of one single 
piece shifts to the ordered behavior of the whole system. There is a huge number of 
phase-locked systems (populations of crickets, yeast cells, lasers ...) and none univer- 
sal mechanism should be expected. The concept of phase pervades the whole physics 
and we will show that its mathematical counterpart touches several intriguing open 
problems. 

Working experimentally it was found that the interleaving of frequencies and 
phases of electronic oscillators interacting in non linear circuits follows arithmeti- 
cal rules. Continued fraction expansions, prime number decompositions and related 
number theoretical concepts were successfully used to account for the experimen- 
tal effects in mixers and phase-locked loops. 1,2 We also made use of these tools 
within the field of quantum optics emphasizing the hidden connection between 
phase-locking and cyclotomy. 3 Finally a class of optimal states in quantum in- 
formation happens to be quantum phase states constructed (phase-locked) from 
Galois fields and rings. 4 In this paper we show that their properties are related to 
incomplete Gauss sums. 

2. Classical Phase-Locking: from Huyghens to the Prime Numbers 

Being obliged to stay in my room for several days and also occupied in making 
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observations on my two newly made clocks, I have noticed a remarkable effect which 
no one could have ever thought of. It is that these two clocks hanging next to one 
another separated by one or two feet keep an agreement so exact that the pendulums 
invariably oscillate together without variation. After admiring this for a while, I 
finally figured out that it occurs through a kind of sympathy: mixing up the swings 
of the pendulums, I have found that within a half hour always return to consonance 
and remain so constantly afterwards as long as I let them go. I then separated them, 
hanging one at the end of the room and the other fifteen feet away, and noticed that 
in a day there was five seconds difference between them. Consequently, their earlier 
agreement must in my opinion have been caused by an imperceptible agitation of the 
air produced by the motion of the pendulums. 

The citation is taken from Ref. 5. The authors cite a later letter by Huyghens 
that the coupling mechanism was in fact a small vibration transmitted through the 
wall, and not movement of air: 

Lord Rayleigh (1907) made similar observations about two driven tuning forks 
coupled by vibrations transmitted through the table on which both forks sat. . . Locking 
in triode circuits was explained by Van der Pol (1927) who included in the equation 
for the triode oscillator an external electromotive force as given in 

d 2 v d 

~dl? ~ ~dl^ 9V ~ P' v3 "> +Uj2v = WoV sinw i, W 

where g is the linear net gain (i.e. the gain in excess of losses, (3 1 the saturation 
coefficient, and u) is the resonance frequency in the absence of dissipation or gain. 
He showed that when an external electromotive force is included, of frequency ujq, 
and tuned close to the oscillator frequency w, the oscillator suddenly jumped to 
the external frequency. It is important to note that the beat note between the two 
frequencies vanishes not because the two frequencies vanish, not because the triode 
stops oscillating, but because it oscillates at the external frequency. 

We can show the locking effect by utilizing the slowly varying amplitude approach, 
including a slowly varying phase $ and oscillation at the external frequency ui a and 
amplitude V 

d$ 

— + if sin $ = oj - luq = u LF , (2) 

where we use u>lf for the detuning term and K — luqVq/V for the locking 
coefficient. 5 

The regime just described is the so-called injection locking regime, also found 
in injection-locked lasers. The equation (2) is the so-called Adler's equation of 
electronics. 6 

One way to synthesize (2) is through to the phase-locked loop of a communica- 
tion receiver. The receiver is designed to compare the information carrying external 
oscillator (RF) to a local oscillator (LO) of about the same high frequency through a 
non linear mixing element. For narrow band demodulation one uses a discriminator 
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of which the role is first to differentiate the signal, that is to convert frequency mod- 
ulation (FM) to amplitude modulation (AM) and second to detect its low frequency 
envelope: this is called baseband filtering. For more general FM demodulation one 
uses a low pass filter instead of the discriminator to remove the high frequency sig- 
nals generated after the mixer. In the closed loop operation a voltage controlled LO 
(or VCO) is used to track the frequency of the RF. Phase modulation is frequently 
used for digital signals because low bit error rates can be obtained despite poor 
signal to noise ratio in comparison to frequency modulation. 7 

Let us consider a type of receiver which consists in a mixer, in the form of a 
balanced Schottky diode bridge and a low pass filter. If /o and / are the frequencies 
of the RF and the LO, and 9{t) and tp(t) their respective phases, the set mixer and 
filter essentially behaves as a phase detector of sensitivity w (in Volts/rad.), that 
is the instantaneous voltage at the output is the sine of the phase difference at the 
inputs 

u(t) = u sm(6(t) - il){t)). (3) 

The nonlinear dynamics of the set-up in the closed loop configuration is well de- 
scribed by introducing the phase difference <I>(i) = 9{t) — ip(t). Using 9 — uj and 
ip(t) = uj + Au(t), with uj = 2nf , u — 2nf and A (in rad. Hz/Volt) as the sensitiv- 
ity of the VCO, one recovers Adler's equation (2) with the open loop gain K = uqA. 
Such a set up is called a phase-locked loop (or PLL). 

Equation (2) is integrable but its solution looks complex. 8 If the frequency shift 
lolf does not exceed the open loop gain K, the average frequency ($) vanishes after 
a finite time and reaches the stable steady state $(oo) = 2ln + sm^ 1 (ujpp/K), I 
integer. In this phase-tracking range of width 2K the RF and the LO oscillators 
are also frequency-locked. Outside the mode-locking zone there is a sech shape beat 
signal of frequency 

u LF = (i>(t)) = (coi F -KY /2 - (4) 

The sech shape signal and the nonlinear dependance on parameters lolf and K 
are actually found in experiments. 8 ' 2 In addition the frequency ujpp is fluctuating 
(see Fig. 1). It can be characterized by the Allan variance ct 2 (t) which is the mean 
squared value of the relative frequency deviation between adjacent samples in the 
time series, averaged over an integration time r. Close to the phase-locked zone the 
Allan deviation is 

*(r) = ^ (5) 

where cto is a residual frequency deviation depending of the quality of input oscilla- 
tors and that of the phase detector. Allan deviation is found independent of r which 
is a signature of a 1// frequency noise of power spectral density S(f) = o-/(21n2/). 
One way to predict the dependence (5) is to use differentiation of (4) with respect 
to the frequency shift ujlf so that 

5Cj lf = 8u; LF (l + ^/CjIf) 1 ' 2 . (6) 
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Fig. 1. Fluctuating counts of the beat frequency (in Hz) close to the phase-locked zone. The 
inputs arc quartz oscillators at 10 MHz. The power spectrum has a pure l/f dependance 



Relation (6) is denned outside the mode-locked zone \u>lf\ > K; close to it, if the 
effective beat note Qlf < K, the square root term is about K/Cjlf- If one identifies 
Sojlf/^lf with a bare Allan deviation er and S&lf/^lf with a magnified Allan 
deviation a one explains the experimental result (5). One can conclude that, either 
the PLL set-up behaves as a microscope of an underlying flicker floor cto, or the 
l/f noise is some dynamical property of the PLL. In the past we looked at a 
possible low dimensional structure of the time series and found a stable embedding 
dimension lower or equal to 4. 9 But at that time the dynamical model of l/f noise 
still remained elusive. 

Adler's model presupposes a fundamental interaction ojlf = \uo — ^(t)\ in the 
mixing of the two input oscillators. But the practical operation of the phase detector 
involves harmonic interactions of the form ljlf — \p^o — Qu(t)\ < lo c = 2irf Cl where 
p and q are integers and f c is the cut-off frequency of the low pass filter. This can be 
rewritten by introducing the frequency ratios v = "^j- and fj, = as \i — q\v— ||. 
This form suggests that the aim of the receiver is to select such pairs (p, q) which 
realize a "good" approximation of the "real" number v. There is a mathematical 
concept which precisely does that: the diophantine approximator. It selects such 
pairs pj_ and qi, coprime to each other, i.e. with greatest common divisor (pi,qi) = 1 
from the continued fraction expansion of v truncated at the index i 



v = {a ;ai,a 2 , ■ ■ -aj = a + l/(di + l/(a 2 + 1/ ■ 



+ 1M)) = -- 

1i 



The diophantine approximation satisfies 

\V < o 



1 



(7) 



(8) 



The fraction ^ is a so-called convergent and the a'^s are called partial quotients. 
The approximation is truncated at the index i just before the partial quotient aj+i. 
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It should be observed that diophantine approximations are different from decimal 
approximations j- for which one gets \v — ^-\ < It was shown in Ref. 1, using 
the filtering condition, that a^+i is given by 



0>i+l 



fo 



fcQi 



(9) 



where [ ] denotes the integer part. For example if one chooses fo = 10 MHz and 
f c = 300 kHz, the fundamental basin ^ = j will be truncated if a i+ i > 33 and 
the basin |± = | w [\\ b e truncated if a i+1 > 6. The resulting full spectrum is a 
superposition of V-shape basins of which the edges are located at 

v\ = {&o; ! 0,1, a,2, ■ ■ ■ , a.i,a,i + i}, 
V2 = {a ;ai,a 2 , ■ ■ ■ ,<k-i, 1, (10) 

where the partial expansion before a,+i corresponds to the two possible continued 
fractions of the rational number 2i. The basin of number v = | = {0; 1,1,2} extends 
to v x = {0; 1,1, 2, 33} = || ~ 0.594, v 2 = {0; 1, 1, 1, 1, 33} = § ~ 0.618. For a 
reference oscillator with fo = 10 MHz this corresponds to a frequency bandwidth 
(0.618 - 0.594). 10 7 MHZ=240 kHz. 

With these arithmetical rules in mind one can now tackle the difficult task of 
accounting for phase-locking of the whole set of harmonics in the beat frequency. 2 

vlf = \p%vo - Qiu(t)\, (11) 

Some essential features can be found in the standard Arnold map model 2 

$„+i = + 2ttO - c sin$„, (12) 

where Q = — is the bare frequency ratio and c = Such a nonlinear map is 
studied by introducing the winding number v = lim„^ 00 ($„ — $o)/(27m). The 
limit exists everywhere as long as c < 1, the curve v versus w is a devil's staircase 
with steps attached to rational values fl = £i and width increasing with the coupling 

Hi 

coefficient c. The phase-locking zones may overlap if c > 1 leading to chaos from 
quasi-periodicity. 10 

The Arnold map is also a relevant model of a short Josephson junction shunted 
by a strong resistance R and driven by a periodic current of frequency cj n and 
amplitude Iq. Steps are found at the driving voltages V r = Rio — r(hujo/2e), r 
a rational number. Fundamental resonances r — n, n integer, have been used to 
achieve a voltage standard of relative uncertainty 10~ 7 . 

To appreciate the impact of harmonics on the coupling coefficient one may ob- 
serve that each harmonic of denominator creates the same noise contribution 
Su>lf = qi5u>(t). They are <j>(qi) of them, where <j>(qi) is the Euler totient function, 
that is the number of integers less or equal to qi and prime to it; the average coupling 
coefficient is thus expected to be l/<f>(qi). In Ref. 2 a more refined model is devel- 
oped based on the properties of prime numbers. It is based on defining a coupling 
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coefficient as c* — cA(n;qi,pi) with A(n;qi,Pi) a generalized Mangoldt function. 2 
It is defined as 

. . . f In 6 if n = b k , b a prime and n = p ; (mod qi), ,„ „. 

A (»;*.ft) = (o otherwise. (13) 

The classical Mangoldt function is A(n) = A(n; 1, 1). It is thus the coupling coeffi- 
cient of the fundamental resonance 1/1. The important result of that analysis is to 
exhibit a fluctuating average coefficient as follows 

1 ' 1 

c£v/c=-£$ZA(n;gi,pi) = ^-y + e(t), (14) 

with e(t) = 0(t~ x l 2 ln 2 (i)) which is known to be a good estimate as long as q t < \ft} 
The average coupling coefficient shows the expected dependance on q t . In addition 
there is an arithmetical noise e(t) with a low frequency dependance of the power 
spectrum reminding 1/f noise. Although that stage of the theory is not the last 
word of the story, it is quite satisfactory that this approach, based on phase-locking 
of the full set of harmonics, is accounting for the main aspects of 1/f noise found 
in experiments. 



3. Quantum Phase-Locking 

Apparently 11 Dirac was the first to attempt a definition of a phase operator by means 
of an operator amplitude and phase decomposition. As we have discussed, with a 
complex c-number a = i?e 4 * one obtains the phase via e 1 * = a/R. Similarly, he 
sought to decompose the annihilation operator a into amplitude and phase compo- 
nents... After a brief calculation we obtain a relation indicating that the number 
operator N and phase operator $ are canonically conjugate 

[JV,<D] = 1. (15) 

The equation immediately leads to a number-phase uncertainty relation which is 
often seen 

SN (5$ > 1/2. (16) 

However, all of the previous development founders upon closer examination. 

This is taken from Ref. 11, a comprehensive review of the quantum phase prob- 
lem. See also Ref. 13. 

To approach the phase-locking problem within quantum mechanics one can start 
from the theory of the harmonic oscillator. The natural objects are the Fock states 
(the photon occupation states) \n) who live in an infinite dimensional Hilbert space. 
They are orthogonal to each other: (n\m) — 8 mn , where 5 mn is the Dirac symbol. 
The states form a complete set: X^^Lo \ n )( n \ = 1- 



a One referee nominated Madelung 12 as the first investigator of quantum phase operators. 
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The annihilation operator removes one photon from the electromagnetic field 

a\n) = y/n\n-l),n= 1,2,- ■■ (17) 

Similarly the creation operator a' adds one photon: a'\n) = \/n + l\n + 1), n — 
0, 1, ■ ■ ■ There is the commutation relation [a, a^] = 1. The operator N = aa 1 * has 
the meaning of the particle number operator and satisfies the eigenvalue equation 
N\n) — n\n). 

Eigenvectors of the annihilation operator are the so-called coherent states \a) 1A 
and are the ones generated by single mode laser operated well above threshold. 

States of well defined phase escaping the inconsistencies of Dirac's formulation 
were build by Susskind and Glogower. 15 They correspond to the eigenvalues of the 
exponential operator 

oo 

E = e^ = (N+l)- 1 ' 2 a= ^\n){n+l\. (18) 

n=0 

Using the Hermitian conjugate operator = e~ 4 *, one gets EE^ = 1, E^E — 
1 — |0)(0|, i.e. the unitarity of E is spoiled by the vacuum-state projector |0)(0|. The 
Susskind-Glogower phase states satisfy the eigenvalue equation E\^f) = e*^|\&); they 
are given as 

oo 

I*) = JV n *|n). ( 19 ) 

n=0 

Like the coherent states the phase states are non orthogonal and they form an 
overcomplete basis which solves the identity operator: j- dtp\ip)(ip\ = 1. The 
operator cos<£> = ^(E + E^) is used in the theory of Cooper pair box with a very 
thin junction when the junction energy E j cos $ is higher than the electrostatic 
energy 16 

Further progress in the definition of phase operator was obtained by Pegg and 
Barnett. 18 The phase states are defined from the discrete Fourier transform (or 
more precisely the quantum Fourier transform since the superposition is on Fock 
states not on real numbers) 

\6 k ) = ^^cxp(2»7r-n)|n). (20) 

V q n=0 q 

The states are eigenstates of the Hermitian phase operator 

9-1 

© 9 = ^0fc|Wfc|, (21) 

k=0 

with 9k = 0q + 2irk/q and #o is a reference angle. It is implicit in the definition 
(20) that the Hilbert space is of finite dimension q. The states \8k) form an or- 
thonormal set and in addition the projector over the subspace of phase states is 
Sfc=ol^fe)(^fel = lg> where 1<? is the unity operator. Given a state \F) one can 
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write a probability distribution |(#fc|-F)| 2 which may be used to compute various 
moments, e.g. expectation values, variances. The key element of the formalism is 
that first the calculations are done in the subspace of dimension q, then the limit 
q — > oo is taken. 18 . 

We are now in position to define a quantum phase-locking operator. 17 Our view- 
point has much to share with the classical phase-locking problem as soon as one 
reinterpret the fraction | in (20) as arising from the resonant interaction between 
two oscillators and the denominator q as a number which defines the resolution 
of the experiment. From now we emphasize such phase states \9' k ) which satisfy 
phase-locking properties and we impose the coprimality condition 

(k,q) = l. (22) 

The quantum phase-locking operator is defined as 

e^ = 5>K><**l. (23) 

k 

with 0k = 27r| and the notation means summation from to q— 1 with (fc, q) = 1. 
Using (20) and (22) in (23) one obtains 

e^ k = ^c,(n-0|n)(Z|, (24) 

n,l 

where the range of values of n, I is from to 0(g), and 4>(q) is the Euler totient 
function. The coefficients in front of the outer products \n)(l\ are the so-called 
Ramanujan sums 

c g (n)=^cxp(2z^-n)^ M y (g) , with qi = q/(q, n). (25) 
k q Wf) 

In the above equation /j,(q) is the Mobius function, which is if the prime number 
decomposition of q contains a square, 1 if q = 1 and (— l) K if q is the product of 
K distinct primes. Ramanujan sums are relative integers which are quasi-periodic 
versus n with quasi-period 4>(q), and aperiodic versus q with a type of variabil- 
ity imposed by the Mobius function. Ramanujan sums have been used for signal 
processing of low frequency noise. 19 In the Ramanujan sum expansion there is a 
modified Mangoldt function b(n) which is the dual of Mobius function 

Kn) = ^A(n) = g^(n,. (26) 

This illustrates that many "interesting" arithmetical functions carry the structure 
of prime numbers. We mention the relation dln((s)/ds = J2 n >i but there 

is also the relation 1/C(s) = J2 n >i ^r- There is a well known formulation 1 of 
Riemann hypothesis from the summatory Mobius function Y^ n =\ ^( n ) = 0(t 1 / 2+e ), 
whatever e. 
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Given a state \(3) one can calculate the expectation value of the quantum phase- 
locking operator as 



<ef k ) = E^I 



(27) 



If one uses the finite form of Susskind-Glogower phase states (19) and a real pa- 
rameter (3 



1 9-1 
— ^exp(m/3)|n>, 

y q n=0 



(28) 



3.5 



2.5 



E 1.5 -i 



! S 




5 10 15 20 25 30 35 40 45 50 
dimension q 

Fig. 2. Oscillations in the expectation value (29) of the locked phase at fi = 1 (dotted line) and 
their squeezing at (i = (plain line). The brokenhearted line which touches the horizontal axis is 
7rA(g)/ In q 



the expectation value of the locked phase becomes 

( ' ock ) = C S - n ) exp(*/3(n - 0)- 



(29) 



For/3 



1 it is found that (6^ ock ) has the more pronounced peaks are at such values 



of q which are powers of a prime number. It can be approximated by the normalized 
Mangoldt function 7rA(q)/lng as shown on Fig. 2. For (3 = the expectation value 
of (Oq° ck ) is much lower. The parameter (3 can be used to minimize the phase 
uncertainty well below the classical value. 3 

Quantum phase-locking effect and its relation to prime number theory has also 
been studied implicitely by Bost and Connes. 20 Instead of an ad-hoc quantum phase 
operator as (21) or (24), it is based on the formulation of a dynamical system 
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and its associated quantum statistics. The dynamical system is first denned by an 
Hamiltonian operator H with eigenvalues equal to the logarithms of integers 

H \n) = lnn|n). (30) 

Using the relations exp(— j3oH n )\n) = cxp(— /3 lnn)|n) = n~^ |n), it follows that 
the partition function of the model at the inverse temperature flo is 

oo 

Trace(exp(-/3 ffo)) - ]T n~ ° = C(/?o), (31) 

n=l 

where C(Po) is the Riemann zeta function 

In quantum statistical mechanics, given an observable Hermitian operator M 
one has the Hamiltonian evolution a t (M) versus time t 

a t (M) = e ltHo Me- ltHo , (32) 

and the expectation value of M is the Gibbs state 

Gibbs(M) = Trace(Mexp(-/?o#o))/Trace(exp(-/3 -ffo)- In Bost and Connes ap- 
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Fig. 3. Phase expectation value (35) in Bost and Connes model at the inverse temperature (i = 3 
(plain line) in comparison to the function /i(q)/c/)(cj) (dotted line) 

proach the observables belong to an algebra of operators \i a and efe which are defined 
by their action on the occupation numbers \n) as 

Ha\n) = \an (mod q)), (33) 
e k \n) =cxp(2HLfc«)|n). (34) 

The first operator fi a acts as a shift in the space of number states; the second one 
ek is such that its action encodes the individuals in the quantum Fourier transform 
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0.01 




-0.07 1 1 1 1 1 1 1 1 1 1 1 

5 10 15 20 25 30 35 40 45 50 

dimension q 



Fig. 4. Phase expectation value (35) in Bost and Connes model at the inverse temperature 
P = 1 + e, e = 0.1 (plain line) in comparison to the function —A(q)e/q (dotted line). 

(20). Like in the quantum phase- locking operator one uses the coprimality condition 
(22) to distinguish in (34) the primitive roots ol unity exp(2iirk / q) , (k,q) = 1. 
One can show that there is a hidden symmetry group which is used to label the 
elements of the algebra b . Using the action of the group, the Gibbs state is replaced 
by the so-called Kubo-Martin-Schwinger (or KMS) state. The system exhibits a 
phase transition with spontaneous symmetry breaking at the inverse temperature 
(3q = 1 which corresponds to the unique pole of the Riemann zeta function C(A))- 
At low temperature /?n > 1 one gets, after tricky calculations, the expectation value 
of the phase operator which replaces (29) in the following form 20 

KMS(eW) = «-A> J] (35) 

p divides q 

p prime 

The KMS state is represented for two limiting cases, the low temperature limit 
ft >• 1 (Fig. 3) and the critical case (3 = 1 + e (Fig. 4), with e ~ 0. In these limits one 
has respectively KMS^»i(g) = ^4 and KMSi+ e ~ ~ A (i) e _ j n the low temperature 
limit the spectrum (26) corresponding to the Ramanujan sum expansion of the 
modified Mangoldt function b(n) = A ^"^^"^ = A(n) is recovered (see (26). Close to 
the critical point (i = 1 + e the oscillations are proportional to A(q) ~ b(q) and are 
of very small amplitude due to the squeezing coefficient e. A comparable squeezing 
effect was already observed in the expectation value (27) of the quantum phase 



b This is the Galois group W = Gal(Q cycl / Q) of the cyclotomic extension on the field of rational 
numbers Q. 
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operator (see Fig. 2). 

Thus after the phenomenological model (14), the Bost and Connes cyclotomic 
model also points to the Mangoldt function as a source of low frequency fluctuations. 
In the last case the model is associated to the spontaneous symmetry breaking and 
the squeezing of phase oscillations at the critical KMS state c . 

4. Galois Phase-locking and Quantum Complementarity 

But what is light really? Is it a wave or a shower of photons? There seems no 
likelihood for forming a consistent description of the phenomena of light by a choice 
of only one of the two languages. It seems as though we must use sometimes the one 
theory and sometimes the other, while at times we may use either. We are faced with 
a new kind of difficulty. We have two contradictory pictures of reality; separately 
neither of them fully explains the phenomena of light, but together they do. " 

(Albert Einstein and Leopold Infold, The Evolution of Physics, p. 262). 

More hints into quantum phase-locking may be discovered by deriving a math- 
ematical view of the complementarity principle. At the conceptual level, two ob- 
servables are complementary if precise knowledge of one of them implies that all 
possible outcomes of measuring the other one are equally probable. The eigenstates 
of such complementary observables are non-orthogonal quantum states, and in any 
attempt to distinguish between them, information gain is only possible at the ex- 
pense of introducing disturbance. 

Mathematically speaking let O be an observable in a Hilbert space of dimension 
q, Tig, which is represented by a Hermitian q x q matrix. Let us assume that its real 
eigenvalues are multiplicity- free and its eigenvectors \b) belong to an orthonormal 
basis B. Let O' be a (prepared) complementary observable with eigenvectors \b') in 
B' . If O is measured, then the probability to find the system in the state \b) E B 
is given by |(6|6')| 2 = 1/q. We here recall that two orthonormal bases B and B' of 
H q are mutually unbiased precisely when |(6|6')| 2 = | for all b e B and b' E B'. It 
can be shown that in order to fully recover the density matrix of a set of identical 
copies of a quantum state, we need at least q + 1 measurements performed on 
complementary observables. As a matter of fact, the mathematical implementation 
of the complementary principle lead us to the search of completes sets of mutually 
unbiased bases (or MUBs for short), a problem which has recently received a peculiar 
attention. 22 

In dimension q = 2, eigenvectors of ordinary Pauli spin matrices (i.e. in dimen- 
sion q = 2) provide the best known example of a complete set of MUBs. It has been 
shown that in dimension q = p m which is the power of a prime p, the complete 
sets of mutually MUBs result from Fourier analysis over a Galois field F q (in odd 
characteristic p) 23 or of a Galois ring R 4 ™ ( when p = 2) (see Refs. 24 and 4 for more 
details). Complete sets of MUBs have an intrinsic geometrical interpretation, and 



c The theory was also used as a model of time perception. 2 
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were for example related to discrete phase spaces 25 or finite projective planes. 26,27 
The complete sets of MUBs in odd power dimension q have a very compact 
Fourier form 4 

\9l) = ^pJ2 ^{ n X an2 + bn )H, a, b e F q . (36) 

in which the coefficient in the computational base {|0), |1), • • • , \q — 1)} represents 
the product of an arbitrary multiplicative character ip(n) by an arbitrary additive 
character n{yn), and where the decomposition y = an + b, has been used. 

Let fix a primitive root g in F q , then every n <G F q is given by n = g s for some 
s€ [0, q — 2] and then 

VW=^-i, (37) 

is a multiplicative character. Every multiplicative character can be obtained in this 
way. 

The additive characters are defined as 

k(x)=J^ x \ ^ p = exp(^), xeF q . (38) 

where tr(x) = x + x p + ■ ■ ■ + x pm 1 is the field theoretical trace. It maps any element 
of F q to an element in the base field F p . 

Eq. (36) defines a set of q bases (with index a) of q vectors (with index b). Using 
a property of Weil sums 4 it is easily shown that, for q odd, the bases are orthogonal 
and mutually unbiased to each other and to the computational base. 

The result of Wootters and Field corresponds to the trivial multiplicative char- 
acter ^o = l- Eq. (36) also defines phase states generalizing those written before in 
(20). The latter are recovered if one uses the trivial additive character n = k in 
(36) d . 



4.1. The phase expectation value 

For the evaluation of the phase properties of a general pure state of an elec- 
tromagnetic field mode in the Galois number field we proceed as in Sect. 3. 
We consider the pure state of the form (28), and we sketch the computation 
of the probability distribution S = | < 9b\f3 > | 2 and phase expectation value 
< ©Gal >= SheF ®b\ < > | 2 • We recall that 9b = 2nb/q (the upper index a for 
the base is implicit and we discard it for simplicity). The probability distribution 
reads explicitly as 4 

— [ ip(—n) cxp(in/3)K(—an 2 — bn)] [ ip(m) exp(— im0)n(am 2 + bra)]. (39) 



d In (20) wc used the (non prime) integer p instead of k. 
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Strictly speaking the phase state defined in (28) is defined only in prime dimension 
q = p. In non-prime dimensions the product in the exponential n/3 of (28) is not in 
F 

For prime dimensions p the phase probability distribution reads 
S = En=i i>( n ) cxp(2i7r(7n + 2an/p)) 
x X)n=i % l ) { m )) exp(— 2i7r(7m + 2am/p)) 
= ^ Z)n,m=i 4>(n)ip{m) exp(2i7r(7(n - m) + a(n - m)(n + m)/p)) 
= E«=i i>(n)i>(n + k) exp(2i7T7fc + a/c(2n + fc)/p) 

= FE£- P+ iexp(2 t 7r7fc)r(fc). 

where we used the notation 7 = —/3/27r + b/p. In the last but one equality above we 
put k = n — m and in the last equality we changed the order of summation, pulling 
out the 7-depcndant factor outside the n summation. The inner sum equals 

min{p— l,p— 1— k} 

T ( k )= ip{n)4){n + k)e-xv(2iirak{2n + k)/p). (40) 

n— max{l,l— k} 

Now, if k = 0(mod p) (which may happen only for k = in the above range) 
the inner sum is trivial and is equal to p. For other values of k the inner sum is at 
most of absolute value 0(p 1 ^ 2 \np) by the Weil bound of incomplete sums 28 ' 29 (note 
that the factor involving 7 is now gone from the sum over n) . 

Hence 

\S(k)\ = \o(l.p + p.p 1 ' 2 \np) = 0{p- x ' 2 hip). (41) 

We get \S\ ~ 0.63 at p=3 and 0.49 at p = 7. Then it decreases slowly with increasing 
dimension p. 

The phase expectation value reads 

<e Ga i>=^ J2 exp(-il3k)T(k)J2 bc M2iKkb/p). (42) 
p k=- P +i 6=1 

The partial sums in the above equation can be evaluated as p(p + l)/2 for k = 
and otherwise 

U = ]T be b = e(l + 2e + 3e 2 + ■ ■ ■ + pe^) = e[-^- - Y~) = ( 43 ) 

where we introduced e = cxp(2iirk/p) and we made use of the relation e p = 1. Easy 
calculations lead to 

\U\ = (44) 

11 2|sin(2fc7r/p)| 1 ; 

An estimate of the phase expectation value can be obtained as follows. 30 Let 
be the smallest (by absolute value) residue of 2fc(mod p). Then | sin(27rfc/p) | = 
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sin(7r|rfc|/p) > 2p\rk\/^ since sin(x) > 2x/n for < x < ir/2. We now define 
rk = 1 for k = 0. Thus, we now have U{k) = 0(p 2 /rk) for any fc. Therefore 

EV- P+1 \u(k)\ = o(p*j:ii- P -^/r k ). 

When k runs between — p — 1 and p+1, with k ^ 0, r k = 2fc(mod p) takes each 
value in the range [0, (p — l)/2] no more than 4 times. The contribution from k = 
is simply 1. So the contribution is 0(p 2 Efet-p-i V^fe) = C(p 2 Ei=i 1)/2 1/*) = 
0(p 2 logp). 

As a result the phase expectation value is bounded by 

I < e Gal > | = ^0(l.p.S^> + p.p 1 / 2 lnp.p 2 lnp) (45) 
= 0(l+p 1 / 2 (lnp) 2 ), 

which is a diverging quantity. 

When V( n ) i s a trivial character equal to 1, the estimate on incomplete Gauss 
sums (40) is replaced by the better bound 0(p/ min (k,p — k)) . It follows that 
S = 0(p- 1 In p)°. 

Here is the proof. When ip(n) is trivial we have T(k) — E^Ljv exp(4i7ra/cn/p) 
where N and M the lower and upper limits in (40). Assume that k ^ 0. Then 
T(k) = (exp(4i7rafciV/p) — exp(4i7rafcM/p))/(l — exp(4iirak/p)). Estimating the 
numerator as 2, we get \T(k)\ < 2/ 1 1 — exp{Aiirak/p)\ = 1/| sin(27rafc/p)|. 

Let now use the notation Sk for the smallest (by absolute value) residue of 
2afc(mod p). Therefore ££-p+i \T(k)\ = O^Efet-p-i V*fc) = 0(plnp) and 
S(A;) = 0(p _1 lnp) as expected. 

4.2. Phase variance 

The phase variance can be written as 

< A9 2 Gal >= £ (0 6 - < 9 Gal >) 2 | < 6 |/ > | 2 . (46) 
fceF, 

For the first term one gets 

V 4 o P+ 1 P 

^0 2 |<0 h |/3>| 2 -— ^ cxp(- i /?fc)T(A:)^6 2 exp(2 i 7rfc6/p). (47) 



6=1 ^ fc=-p+l 6=1 



3 



The partial sums in the above equation can be evaluated as ^- for k = and 
otherwise 

Using the same type of reasoning than in the last section 

*£0(l.p.£ + p.p 1/2 lnp.p 3 lnp = 0(1 + 3p 1 / 2 (lnp) 2 ). (49) 



e However one way to squeeze the bound on the phase expectation value to ir is to project on the 
state |/3) = |0) in (28) as it was done in Fig. 2 in the context of "Ramanujan-type" phase states . 
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The second term is 



p , p+i p 

< e Ga i > 2 J2\ < °b\P > \ 2 = — eM-iPk)T(k) exp(2inkb/p). (50) 

6=1 k=-p+l 6=1 

The inner sum equals p if k — and otherwise, so that the whole contribution os 
0(1). 

Finally the third term is 

p 

-2 < e Gal > J2 °b\ < e b \P > | 2 - -2 < e Gal > 2 . (5i) 

6=1 

Using (46) the absolute value is is bounded by 

0(l+p{\npf). (52) 
All estimates of the contributing terms in the variance are diverging with p. 

5. Maximally entangled states 

By definition entangled states in Tt q cannot be factored into tensorial products of 
states in Hilbert spaces of lower dimensions. There is an intrinsic relation between 
MUBs and maximal entanglement. 

Generalized Bell states may be defined using the multiplicative Fourier transform 
(21) applied to the tensorial products of two qudits, 

1 9-1 

|B u>fc ) = — J2uj k q n \n,n + u). (53) 

V q n=0 

Also these states are both orthonormal, {B u> k\B u ' ,k') — ^uu'^kk'i an d maximally 
entangled, trace2\B u .k){B u ^k\ = \lq- A more general class of maximally entangled 
states is obtained using the Fourier expansion (36) over F q (q odd) as follows 

1 ^ 

\Kb) = -y= E 4 rlian+ " )n] l n ' » + «> ■ (54) 



71=0 

2 



In general, for q a power of a prime, starting from (54) one obtains q bases of q 
maximally entangled states. Each set of the q bases (with u fixed) has the property 
of mutual unbiasedness. Similarly, for sets of maximally entangled m-qubits one 
uses the Fourier transform over Galois rings. 4 



6. Conclusion 

The phase relation between a single piece and the whole system is strongly contex- 
tual, but the working mathematics is amazingly universal. In an electronic phase- 
locked loop we found that the position of mode-locked zones is controlled by the 
arithmetic of irreducible fractions (7)-(ll), and the strength of lockings is related to 
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prime number theory via the Mangoldt function (14). We also developed two differ- 
ent approaches of phase-locking within the context of quantum optics. One of them 
uses a discrete phase operator (23), and the phase expectation value also relates 
to prime number theory via Ramanujan sums (see (27)). In a more sophisticated 
form it is linked to the Riemann zeta function (see (35)). In a second approach, the 
phase locked states are properly defined Fourier transforms over a Galois field (see 
(36)). They connect to mutually unbiased bases, which appears in the mathematical 
formulation of quantum complementarity Incomplete Gauss sums are at the kernel 
of phase variability in this case. 
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